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Abstract: This paper presents a novel controller design for half-car magneto-rheological (MR) 
suspension systems. By introducing a piecewise approximation model, the nonlinear constraints 
of MR dampers are transformed into piecewise constant constraints. Further, by defining a 
piecewise affine control law, the attenuation of pitch and heave responses are realized by 
satisfying a given H^ index. The control gains are derived through linear-matrix-inequality 
(LMI) optimization. Simulations under random and impact roads have verified the effectiveness 
of proposed controller on performance improvement and constraint satisfaction. 
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1. INTRODUCTION 

A typical automotive suspension must provides ride com¬ 
fort in view of road roughness inputs and secures good 
road handling property. Conventional passive suspensions 
always suffer conflict from these two requirements, which 
has resulted in massive researches and attempts in the area 
of active and semi-active suspensions (Tseng and Hrovat 
(2015)). Active suspensions can offer opportunities for 
substantial improvements in ride and handling as well as 
overall vehicle posture, but high energy consumption, cost 
and complexity are requested with pumps or compressors 
(Huang et al. (2018)). Semi-active suspensions can improve 
the comfort and handling performance by only adjust 
the damping properties with acceptable energy consump¬ 
tion and compact actuator size, which are applicable and 
promising. 

Among the existing semi-active actuators, magnetorhe- 
ological (MR) dampers, which can change its damping 
with different applied current, have drawn more attentions 
from the academic and industrial circles due to their fast 
response and highly adjustable range (Choi et al. (2016)). 
However, many researches have shown that MR dampers 
have complex dynamics, for example, the hysteresis, bi- 
viscous and energy dissipativity (Wang and Liao (2011)). 
Thus, the task of controller design while considering the 
actuator dynamics has been a hot issue. 

Skyhook control or its modified styles have been widely 
used in the industrial area (Choi et al. (2001)), where 
an ideal damper that connects the vehicle body and the 
‘sky’ is set as a desired closed-loop property. To further 
explore a control method with higher performance, many 
control theories, such as LQR/LQG, H oo, sliding mode 
and so on, have been employed. For example, Du et al. 
(2005) has firstly design a static Hqq controller for an 
active suspension, and then the active control force is 
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realized with the MR damper using a designed polyno¬ 
mial model. In Pohoryles and Duffour (2015), the active 
damping force is designed by LQG control, and a new 
rule to obtain the current value is derived by analyzing 
the dynamics of MR dampers under different excitation 
frequencies. However, these methods have neglected the 
dynamics of MR dampers during the design of damping 
force, thus the proposed results are ‘clipped optimal’. To 
consider the complex actuator nonlinearities, some control 
algorithms have been given based on the MPC technique. 
For example, by designing a hybrid MPC controller to ap¬ 
proach the optimal semi-active suspension, Giorgetti et al. 
(2006) has shown that significant deviation of optimal 
from the ‘clipped optimal’ can occur. To realize the online 
application, Canale et al. (2006) and Csek et al. (2015) 
have explored a ‘fast’ MPC method and an neural network 
approximation of MPC laws, respectively. 

Considering the above problems, a static H^ controller 
which can be easily applied has been given in this paper 
based on a piecewise approximation of the suspension 
model. Additionally, differ from the above MPC meth¬ 
ods which only consider a quarter-car model, this paper 
studies the attenuation of pitch and heave motions with 
a half-car model. The nonlinear actuator constraints are 
approximated by piecewise constant constraints and the 
system state space is divided by a certain number of par¬ 
titions based on the vibration velocities of front and rear 
suspensions. This paper is organized as follows. Section 2 
gives the half-car suspension model and the nonlinear actu¬ 
ator constraints. The model approximation and controller 
design are given in Section 3. Section 4 gives the simulation 
results and analysis under different road profiles. 

2. PROBLEM STATEMENT 
2.1 Half- Car Suspension System 

A half-car MR suspension model is shown in Fig. 1 where 
Ffmr and F rmr are damping forces produced by the MR 
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Fig. 1. Half-car MR suspension systems 

dampers. The dynamics of the sprung mass (M s ) and 
unsprung masses (m f u and m ru ) can be described as follows 

M S Z S = .Ff mr F rmr Ff s F rs Ff ( \ F r( \, 

h<P = {F fd + Ffe + Ffmr) — ^r(F" r d + F rs + T rmr ), q\ 

^fulfil — Ff mr T Ffs T Ffd Fft, 

^ru^ru = -F rmr + F rs + F rc f Frt 

where z s is the vertical displacements of vehicle mass 
center, I s and ip stand for the moment of inertia and the 
pitch angle of vehicle body. The forces in equation (1) 
produced by springs, passive dampers and tires satisfy 

Ffs = ^fs(^s If ^fu)5 F'rs = k rs (z s + l Y (f Zru)i 
Ffd — Cf (is IfP ifu)j Frd — ^r(is T ^r^ iru)} (2) 
Fft — ^ft(^fu ^fr)? Frt — ^rt(^ru ^rr) 

where fcf s and k rs are the spring stiffness; Cf and c r are 
the damping coefficients of passive dampers; kf t and k rt 
are the equivalent stiffness of pneumatic tires; Zf u and z ru 
stand for the vertical displacements of unsprung masses; 
Zf r and z rr are the excitations from uneven road profiles. 


2.2 Nonlinear Actuator Constraints 

In equation (1), Ff mr and F rmr have complicate relation¬ 
ship with the driving current I and the deflection condi¬ 
tions of MR dampers, which can be generally described by 
a modified Bouc-Wen model as follows 

Fmr = T ki(x — Xo)? 

i = ———[azu + coX + fco(x-0]> ( 3 ) 

c 0 + Cf 

™ = ~i\x - £IM n_lcc7 - P(x - OM" + Mx - 0 

where F mr and x are the damping force and stroke of 
MR dampers, respectively (referring to (1), F mr = Ff mr 
or F rmr and x — IfP %fu or z s T l x p z rn )i C 

and w are the inertial variables. Furthermore, model (3) 
has ten parameters which form a parameter set <£> = 
[co, ci, &o, fci, k 2 , <a, /3, 7 , xo? n \ where <a, Co and c\ satisfy 
linear or nonlinear relationships with 1. 

Due to the complexity shown in (3), it is difficult to design 
a controller directly based on a model consists of ( 1 ) and 
(3). Therefore, a practical method generally consists of 
two steps. In the first step, a controller is designed based 
on ( 1 ) to give desired damping forces, and in the second 
step, the forces are transformed into current values based 
on calculations of MR damper models or static MAPs. 
However, in the first step, the capacity of MR dampers 
should be totally concerned. Otherwise, the desired forces 
will be impossible to be tracked through the second step, 
and the designed performance will not be guaranteed. 



Fig. 2. Force/velocity characteristics of MR dampers 

Fig. 2 shows the simulation results of model (3) under 
sinusoidal excitations and varied current values, where the 
nonlinear dynamics of MR dampers, including hysteresis, 
bilinear and energy dissipation characters, are revealed 
under the force/velocity coordinate system. Meanwhile, 
Fig. 2 also shows the available scope of the damping 
force. In this paper, the envelope of the available scope 
is described with six straight lines as follows 

9k(x) = a k + b k x, k = 1, (4) 

Further, by defining the endpoints of these lines as v^i = 
1,..., 4, the bounds of the available scope can be obtained 
as follows 

{ 9i (or 4 ) (x) Vl < X < V 2 

52 (or 5) (x) V 2 < x < V3 (5) 

93 (or 6) (x) V 3 <x<v 4 

where g u stands for the upper bound and g\ is the lower 
bound. Based on (5), the nonlinear constraints on Ff mr 
and F rmr can be described as follows 

9 \(zs IfP ifu) 5s Ff mr 5 0u(is IfP ifu)5 /0\ 

01 (is F l x p iru) 5 ^ F rmr ^ 0u(is T l x p iru)* 

2.3 Problem Statement 


Define the state vector of system as 

X = [Xi X 2 Xs X/f £5 X 6 X 7 X 8 ] T , 

X\ — Z s If p £f u , X2 — Z s T l r p £ m , 

X3 — Zf u Zf r , X4 — Z ru Z rr , X§ — Z s If Pi 

X 6 =Z s + l r (p X 7 = Zf u , Xg = iru- 


Then, by defining the control input as u = [Ff mr F rmr ] T , 
the disturbance input as w = [if r i rr ] T and the perfor¬ 
mance output as z = [pi'z s p 2 p] T 1 expression (1) can be 
written in a state space equation as follows 
x = Ax + Bu + Ew , 

2 : = Cx + Du ^ ' 

where 

A = 
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Fig. 3. Piecewise approximation of the available scope 
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Further, the nonlinear constraints in (6) can be written as 
51 ( 2:5 - x 7 ) <ui< g u (x 5 - x 7 ), 

51 ( 2:6 - a: 8 ) < u 2 < g u (x 6 - a: 8 ). 


(9) 


Above all, the main task of this paper is to design control 
inputs Ui (i = 1 , 2 ) such that system ( 8 ) is asymptotically 
stable and the H^ norm from w to z is less than a 
given scalar 7 . In addition, the nonlinear state-dependent 
constraints in (9) should be satisfied. 


3. PIECEWISE AFFINE CONTROLLER 
SYNTHESIS 


To handle the nonlinear constraints in (9), existing re¬ 
search always appeal to MPC techniques which face large 
calculation burden. Therefore, this paper will consider a 
static iLoo control method based on a piecewise model 
approximation technique. 

3.1 Transformation to Constrained Piecewise Systems 

As shown in Fig. 3 which takes the MR damper of front 
axle as a example, the available scope can be approximated 
by 9 strip-shaped partitions (3 partitions for each of the 
three different damping regions: pre-yield region, post¬ 
yield regions 1 and 2). Define i = 1,..., 9 as the partition 
index of front axle, pm and pui as the lower and upper 
velocity bounds of the i-th partition, then we have 


1 

f Vi 

+ 

(v 2 - 

- Vi )(i - 

l)/3 

1 

< i 

< 3 


PLi = < 


+ 

(V3 - 

- v 2 )(i - 

4)/3 

3 

< i 

< 6 

(10) 

1 

1 

l V 3 

+ 

(v 4 - 

-v 3 )(i- 

7)/3 

6 

< i 

<9, 


1 

f Vi 

+ 

(V2 - 

-vi)i/2, 


1 

< i 

< 3 


PUi = \ 

1 ^2 

+ 

(V3 - 

- v 2 )(i - 

3)/3 

3 

< i 

< 6 

(11) 

\ 

l V3 

+ 

(v 4 - 

- v 3 )(i - 

6)/3 

6 

< i 

< 9. 



Fig. 4. System partitions and state transition 

Further, we use two constant values to approximate the 
nonlinear control constraint in each partition. Specifically, 
by define qm and qm as the damping force constraints of 
i-th partition and referring to (5), (10) and (11), we can 
obtain the expression for qm and qm as follows 


9\{PLi) 

i ± 5 

(12) 

(g\{PLi) - 5u(P!7i))/2 

i = 5, 

9u(pui) 

i + 5 

(13) 

(. 9u(PUi ) - 9\(PLi))/2 

i = 5. 


Based on the above approximation, the nonlinear con¬ 
trol constraint of front axle has been transformed into 
piecewise linear constraints. Furthermore, benefiting from 
the symmetry of the actuator dynamics in the first and 
third quadrants of the velocity/force coordinate system, 
we can choose the values of v 2 and U 3 such that v 2 = — U 3 . 
Then, we can guarantee that the system origin (x = 0) 
will not located in any boundaries of the partitions and 
the constraints in origin-included partition {i = 5) are 
symmetrical, thus qm = 0U6- 

Similarly, by defining j = 1,..., 9 as the partition index of 
rear axle, we can obtain the velocity bounds ppj and puj , 
and the corresponding control constraints qLj and qup 

By using the above approximation technique and choosing 
appropriate sample time and a discretization method, ( 8 ) 
and (9) can be transformed into a combination of the 
following discrete-time subsystems 

x(k + 1 ) = Aijx(k ) + Biju(k ) + Eijw(k ), 

z{k) = Cijx(k ) + Diju(k ), (14) 

u(k) G Sij , x(k) G 7 Zij 

where IZij := {x(k)\p Li < x&(k) - x 7 (k) < PufVl) < 
xq( k) — xs(k) < Puj} is the defined state domain and 
Sij := {u(k)\q Li < ui(k) < qui.QLj < u 2 (k) < quj} is the 
control constraint of the (i,j)-th subsystem. 

Fig. 4 shows the state domain of all subsystems. It is easy 
to see that there will be 81 subsystems when the available 
scope of each MR damper is divided into 9 partitions. 
Further, Fig. 4 has also taken the (i, j)-th domain as an 
example to show the possible state transition. Thus, the 
boundary of any two domains where the state transition 
may happen is assumed to be a line (L-bounds in Fig. 4) 
but not a single point (P-bounds in Fig. 4). In applications, 
when the state unfortunately pass through the P-bounds, 
this assumption is easy to be realized by adding a small 
offset to the state value. In addition, when the system 
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state transits from 7 Zij to domain 7^ m?n at the time k , we 
specify that the system is governed by the dynamics of the 
subsystem in 7 Zij at that time. Above all, let 17 represent 
the set of index pairs which denote the possible switch 
from one partition to itself or another partition, that is 

n = ( m,n))\x(k ) G 1lij,x(k + 1) G Hm,n\ (15) 

Therefore, the control problem in section 2.3 has been 
transformed into: design control inputs Ui(i = 1,2) such 
that subsystems in (14) are stable and the H^ perfor¬ 
mance is satisfied for all of the state transition defined in 
(15). In addition, the constant constraints in each partition 
should be satisfied. 

Remark 1. The approximation quality can be enhanced 
by employing more partitions, however the system and 
control design will become more complex. So, this paper 
has taken 9 partitions for each MR damper as an example, 
but the approximation technique and the following control 
method are not limited to this. 

3.2 Design of a Piecewise Affine Controller 

From Fig. 3, it is obviously that the control constraints 
in (14) are mostly asymmetric about zero (except for the 
origin-included partition 7 ^ 5 ), thus qLi —Qui and qLj 
quj. However, asymmetrical constraints are always hard 
to be handled since most of the synthesis methods rely 
on the concept of invariant set which are ellipsoidal and 
symmetric. Consider this problem, we employ a piecewise 
affine control law as follows 

u(k) = K itj x(k ) + ij = 1,...,9 (16) 

where the affine items are pre-designed as = [(qu + 
qui )/2 (qLj + quj)/2 ] T • Substitute (16) into (14), we can 

get the closed-loop system as 

x(k + 1) = Ai nx(k) + Ei nw(k) + di 7 -, 

z{k) = Cijx(k) + eij, (17) 

I ^D-ijxfkf) ) (<r) I — x(k) G IZij 

where Aij = A{j T Fid R id , G id = G i,j A Dij R ij , 
dij = BijCij , eij = DijCij, a G {1, 2} is the index of row 
elements, u iJ(y 1 ) = (q m -qLi)/% and Uij( 2 ) = ( Quj~qLj )/2. 
Therefore, the asymmetrical control constraints has been 
transformed into the symmetrical ones on Kijx. Further, 
as constant items of output equations have no influence 
on the system stability and FL^ performance, the control 
gains Kij can be equally designed based on the following 
system 

x(k + 1) -• A it jx(k) + Eijw(k ) + dij, , , 
z(k) = Cijx(k). 

The definition domains 7 Zij can be expressed by ellipsoids 
as follows 

Cij = {x\\\FijX + fij\\ <1}, i,j = (19) 

where F iji(1) = 2M {1) /(p m - p Li ), F i>j(2) = 2M^)/(puj ~ 

PLj ), fij( 1 ) = -(PLi +PUi)/(PUi ~ PLi ), fij( 2 ) = -(PLj + 
puj)/{puj - PLj), Fij { 1 ) (or fij ( 1 )) and F i)j(2) (or 
are the first and second row elements of Fij (or fj ), 
respectively, and 

M (1) = [0 0 0 0 1 0 - 1 0], 

M (2) = [0 0 0 0 0 1 0 - 1]. 


By specifying the energy bound of w(k) as a known 
constant, thus ££T 0 ||ie(&) || 2 = w m ax, a theorem to obtain 
Kij can be derived as following. 

Theorem 1 . For the given scalars 7, re max and all the index 
pairs ((i, j), (m, n)) G 17 , system ( 18 ) is asymptotically 
stable, the constraints on Kijx are satisfied and the 
norm from w to z is less than 7, if there exist matrices 
QiJ ^ Qm,n Qm,n ^ G id ' 1 R id and SCalarS 

rjij > 0, fiij > 0 satisfying the following inequalities 
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d L * 
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Ti fJ - 


5 , ( 20 ) 
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E id 
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L i,3 

Y 

^ i,3 -1 


( 21 ) 


b,j(cr) 


7 7C max 

V^(cr) R id 

F G 

r '^,3 Ksr ^,3 


Fj * 
-I 


★ 

* 


0 fij 1). 


r 7/ 2 . 

9 1,1 * 

7 TC max 

L —I] 


< 0, i,j = 5 


<0,i,j/5, (22) 


(23) 


where ^ij — {AijGij + BijRij ) T , T^j — {CijGij + 

D ij R ij) T i = Qid - G iJ - G h' ^(1) = I 1 °] and 

^( 2 ) = [0 1]. Then, the control gains can be derived as 

Kid= R id G i]- (24) 


Proof. Due to the limited length and space, the proof 
is omitted here. Similar techniques used to conclude this 
theorem can be referred in Gao et al. (2009) where a 
77oo controller is designed for piecewise affine systems with 
symmetrical control constraints. 

Remark 2. In theorem 1, (20) and (21) are used to guar¬ 
antee the Hqq performance and system stability, while 
(22) and (23) are designed for the satisfaction of control 
constraints. Further, as (20) and (22) can only be used 
when fijfJj — I > 0 which is only satisfied for origin- 
excluded partitions, (21) and (23) have to be specially 
designed for the origin-included partition. 

Remark 3. Based on theorem 1, for the system partitions 
defined in section 3.1 and the state transition shown in 
Fig.4, there are totaly 369 inequalities to be solved. This 
task can be handled by employing a Matlab toolbox named 
SeDuMi (Sturm (1999)) together with a optimization 
software named MOSEK. Finally, according to (16), we 
will obtained an independent affine controller for each of 
the 81 partitions. 
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heave impact pitch impact 



random excitation wavelength < Z f + l r I random excitation j wavelength > 2(1 f +l r ) random excitation 


Fig. 5. Simulation road profile consists of random excitations, heave and pitch impacts. 


4. SIMULATION ANALYSIS 
4-1 Parameters and Environments 

To verify that the proposed piecewise affine Hqq (PWA- 
iLoo) control design achieves the desired objectives, the 
resulting closed-loop system has been simulated under 
different road conditions and compared to both a passive 
suspension system and an semi-active suspension con¬ 
trolled by classic Skyhook controller. In our simulations, 
the parameters of the half-car suspension model are listed 
in Table 1. 

Table 1. Parameters of half-car model 


Parameter 

Value 

Parameter 

Value 


575kg 

h 

1.38m 

m fu = m ru 

60kg 

lr 

1.36m 

&fs “ ^rs 

16812N/m 

Is 

769kg/m 2 

Cf = C r 

500N/m 

k ft — krt 

190000N/m 



(a) 


(b) 


Fig. 6 . PSD results under random road profiles 


Further, to describe the nonlinearities of MR dampers, a 
hyperbolic tangent function-based model from Guo et al. 
(2006) is employed. The model can be expressed as follows 

F mr = C t ox + ktox + Vmv I tanh(c t ix + hix)- (25) 

According to the test results of a real MR damper from 
FAW-Tokico Corp., the parameters in (25) can be ob¬ 
tained as Cto=810.78, k t o=620.79, c t i=5.78, k t i=20, and 
2/mr=1200. In the simulation, we assume that Ff mr and 
F rmr both conform to (25). Then, by choosing —v\ = 
t’ 4 = 3 m/s and —V 2 = t’ 3 = 0 . 4 m/s, we can obtain the pa¬ 
rameters of the envelope lines in (4) as follows 

a\ = 42.31, a 2 = 1836, <23 = 3635, a± = —3508, 

<25 = —1762, <26 = —5.989, bi = 825, 62 = 5240, 

b 3 = 800.6, 64 = 859.8, b 5 = 5202, b 6 = 827.3. 


The desired damping forces (Ff mr and F rmr ) of PWA-H^ 
control are obtained according to (16) and for Skyhook 
control, the desired damping forces are given as 

Lfmr — Cfsky^s ^f^)j -Grmr — ^rsky^s T Iv^P) (26) 
where Cf s k y = c rs k y = 4000Nm/s. Both the damping forces 
of PWA-iLoo and Skyhook controllers can be transformed 
into the current values with the inverse expression of (25). 
Details of this transformation are omitted in this paper. 


A excitation designed for the simulation is shown ir 
Fig. 5 which consists of random excitations and impacl 
excitations. The expression for this excitation is as follows 

l G [0, 20) 
l G [20, 20 + L) 
l G [20 T L, 40 T L) 
l G [40 T T, 40 T 3T) 
l G [40 + 3L, 60 + 3 L 


^fr(0 — < 


ISO c 

A m i sin((/ — 20)/L) 
ISO c 

dm 2 sin((Z - L - 40)/2 L) 
ISO c 


where L = Zf + Z r , l is the longitudinal displacement, l = vt , 
v is the vehicle speed and A m i (i = 1 , 2 ) are the amplitude 
of impact roads. For the first impact road, we choose 



Fig. 7. Control output and state trajectory 

the wavelength equal to L such that the heave motion 
will be the major response. For the second impact road, 
the wavelength is set as 2 L to trigger a pitch response. 
The random road profiles will trigger both heave and 
pitch motions and are chosen as the class C road which 
is suggested by the international standard organization 
(ISO) and are simulated by the filtering of a white noise. 
In this paper, we choose v = 50km/h and A m i = A m2 = 
0 . 2 m. 


4-2 Results and Analysis 

Fig. 6 shows the estimated PSD results of the heave and 
pitch responses under random road. It has been clearly 
shown that the proposed PWA-iLoo controller has evident 
advantages in a frequency range of l-4Hz when compared 
with the Skyhook controller. However, in the higher fre¬ 
quency band, the performance of PWA-iLoo controller is 
poorer than that of the Skyhook controller, but improve¬ 
ment with respect to the passive one can still be achieved. 
Fig. 7 shows the control outputs of different controllers and 
the transfer of system state under PWA-LGo control. Refer 
to Fig. 7(a), we can see that the output of the proposed 
PWA-iLoo controller can satisfy the nonlinear actuator 
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Fig. 8. Responses under heave impact. 



Fig. 9. Responses under pitch impact. 



Fig. 10. Control output ((a):heave impact, (b) pitch im¬ 
pact ) 

constraint, while that of the Skyhook controller are mostly 
unavailable. In Fig. 7(b), the system states under PWA- 
Hqq control transfer among the designed system regions. 
However, a insight into the dense regions shows that the 
state switches frequently occur in the region boundaries, 
which may influence the system performance and need 
further analysis in the future works. 

Figs. 8 and 9 give the responses of different controllers 
under heave and pitch impact roads, respectively. Different 
from the random road simulation where the ie max is only 
chosen as 20, we simulate the conditions when re max is 
set as 20 and 200 for the impact roads. From Figs. 8 
and 9, we can see that the proposed PWA-LGo controllers 
have smaller peak values and shorter settling time than 
the Skyhook controller and the passive suspension. Ad¬ 
ditionally, the performance of PWA-LGo controller when 
tiimax = 200 is slightly better than that when re max = 20. 
This phenomenon can be explained by referring to Fig. 10 
where the control output of the two PWA-LGo controllers 
are given and the output of PWA-LGo when re max = 20 
have located beyond the available region of MR dampers. 
Above all, we can conclude that although the nonlinear 
constraint can be handled using the proposed method, but 
using a single controller to adapt different road profiles will 
also result in conservatism. 


5. CONCLUSION 

A novel method to handle the nonlinear actuator of 
half-car MR suspensions has been given in this paper. 
The nonlinear control constraints have been approximated 
by piecewise constant constraints and the state transfer 
among different system partitions has been defined. Based 
on the new modeling method, a PWA-LGo controller 
has been given. The simulations under random roads 
and impact roads have revealed the effectiveness of the 
proposed method. Future works will focus on the smooth 
transfer of system state between different partitions and 
the gain scheduling method for varied road conditions. 
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